A multi-blob representation of semi-dilute polymer solutions 
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A coarse-grained multi-blob description of polymer solutions is presented, based on soft, transfer- 
able effective interactions between bonded and non-bonded blobs. The number of blobs is chosen 
such that the blob density does not exceed their overlap threshold, allowing polymer concentra- 
tions to be explored deep into the semi-dilute regime. This quantitative multi-blob description is 
shown to preserve known scaling laws of polymer solutions and provides accurate estimates of am- 
plitudes, while leading to orders of magnitude increase of simulation efficiency and allowing analytic 
calculations of structural and thermodynamic properties. 
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Many conformational, structural and thermodynamic 
properties of semi-dilute polymer solutions, both in the 
bulk and under confinement, can be qualitatively un- 
derstood in terms of scaling arguments based on the de 
Gennes-Pincus "blob" picture . This picture is of value 
whenever the characteristic length scale of the polymer 
solution (e.g. the correlation length ^), or of the con- 
finement is significantly shorter than the radius of gyra- 
tion Rg of an isolated polymer chain. The blob picture 
suggests a systematic coarse-graining procedure whereby 
each polymer chain is divided into a number n of blobs, 
each containing the same number of monomers or seg- 
ments of the initial coil, such that blobs of the same or 
different chains do not, on average, overlap. In this let- 
ter we present and validate a quantitative formulation of 
such a multi-blob representation, which allows a popular 
single blob coarse graining procedure 0, 0, 0, H, Q to be 
extended to highly concentrated solutions. 

Consider a solution of N self avoiding polymer chains 
of L segments (each of size b) in a volume V; the polymer 
number density is p = N/V, and if Rg bL'^ (with 
ly ~ 0.588 the Flory exponent) is the radius of gyration 
of an isolated chain, the polymer overlap density is p* — 
3/ (47ri?3) ^ 

In the dilute regime, p* < p, where polymers do not, 
on average, overlap, they may be represented by a single 
blob of radius Rg] the effective interaction potential v{r) 
between the centres of mass (CM) of two blobs can be 
calculated by averaging over monomer conformations for 
a given distance r between their CM's, e.g. by Monte 
Carlo (MC) simulations of an isolated pair of polymers 
0, 0, [j, H, The resulting v{r) depends weakly on 
polymer length L, and in the scaling limit (L — > oo), it is 
accurately represented by a Gaussian of width ~ iig : 



v(r) 



Aexp [-air/Rg)"^ 



(1) 



where A ~ 1.80 and a ~ 0.80; the softness of the re- 
pulsive interaction, characterised by a modest free en- 
ergy penalty of ~ 2kBT at full overlap (r = 0) of two 



polymers, reflects the low average monomer concentra- 
tion c ~ L^^'^'^ ^ L^"-^^ inside each coil for long chains. 

In the semi-dilute regime, p > p* , polymer coils 
overlap, and this is reflected in a significant density- 
dependence of the effective interaction (H, which spoils 
the simplicity of the coarse-graining procedure and in- 
troduces complications associated with state-dependent 
interactions 7]. This density-dependence signals the 
fact that in the semi-dilute regime the relevant length 
scale is no longer Rg, but the shorter correlation length 
i - Rg{p/p*)-^, with 7 = - 1) - 0.77 [H. 

These shortcomings may be overcome by switching to 
a multi-blob representation, where each of the n blobs 
is made up of / = L/n segments. If rg ^ bl'^ is the 
blob radius of gyration, the blob overlap concentration 
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pl = 3/ (47rr3) = p*, 
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This means that the poly- 



mer density p — pb/n can be increased beyond p* , up 
to n^'''~^p* ~ rt° ''^p* before the blobs overlap. In other 
words the more blobs are used to represent one polymer, 
the deeper one can penetrate into the semi-dilute regime 
without significant blob overlap. Under those conditions 
the effective interactions between the CM's of the blobs 
are expected to be practically independent of blob den- 
sity ph, and may be taken equal to their zero density 
limit. 

These effective interactions include the pair potential 
v{r) between non-bonded blobs, and the "tethering" po- 
tential (/)(r) between adjacent blobs on a given chain. The 
former is expected to be similar to the Gaussian repul- 
sion in equation ([T]), with Rg replaced by rg, i.e. the 
same as the effective potential between polymers in a 
single blob representation [1, 0, @]. 0(r), on the other 
hand, may be expected to be the superposition of v(r) at 
short distances r, and a harmonic spring at large elonga- 
tion, similar to the entropic spring of a Gaussian chain 
[l|,[g|, albeit with a renormalized spring constant: 



(r) = v{r) + -(r - ro) + c. 



(2) 



These conjectures are borne out by MC simulations of an 



2 




FIG. 1: (Color online) Effective potentials (l>{r) between 
bonded (right hand scale, continuous line) and v(r) between 
non-bonded (left-hand scale, crosses) blobs in units of ksT, as 
computed by MC simulations of an L = 2000 SAW polymer, 
divided into 2 equal blobs. v{r) is accurately fitted by the 
Gaussian ((T} with A = 1.752 and a = 0.804. 0(r) is repre- 
sented by the form ((2)1 (dashed curve) with the harmonic teth- 
ering potential fitted by 0.534(r/rg - 0.731)^ - 0.631 (dots). 



isolated pair of self-avoiding walk (SAW) polymers on a 
cubic lattice; each polymer is divided into two sub-chains 
of equal length L/2, and the intra-molecular distribution 
function s{r) of the CM's of the two blobs, as well as 
the inter-molecular pair distribution function g{r) of the 
CM's of blobs on different chains are monitored as func- 
tions of the CM-CM distances r. (f){r) follows directly 
from s(r) ^ exp(— 0(r)/fcBr), while v{r) is related to 
g{r) and s{r) via an exact inversion relation [9|, which 
we used previously for simple models of a di-block copoly- 
mer fiol . The resulting 0(r) and v{r) are shown in 
figure 1. v{r) is indeed of Gaussian form ([T]), and is prac- 
tically indistinguishable from the low density limit of the 
effective pair potential between the CM's of two polymers 
of the same length, within the single-blob case [5|, [y]. 

So far no approximation has been made, and the effec- 
tive potentials in figure[T]are "exact" , for a given length L 
of the polymer. As in the single blob representation [aQ, 
the L-dependence of v{r) and 0(r) is found to be negli- 
gible for L > 1000. We now make the reasonable trans- 
ferability assumption that the interaction potentials in 
figure [H as functions of the scaled distance r/rg, remain 
identical, whatever the number n of blobs into which a 
polymer chain is divided, provided the number / of seg- 
ments in each blob remains sufficiently large. Note that 
within this assumption, the potential ^ acts between 
adjacent (tethered) blobs only, while v{r) is the same for 
all pairs of blobs of the same, or different, chains. This 
defines the multi-blob model of homopolymer solutions, 
and applies as long as pb < pI- Note that in view of the 
"softness" of the effective interactions the model is par- 
ticularly well suited for simulation purposes (apart from 
the obvious speed-up due to the considerable reduction 



of the number of degrees of freedom from 3N-L to 37V • n) 
as well as for approximate theories 13, to which we 
shall return. 

We now put the above multi-blob representation to the 
test, focussing on the structure factor S{k) of a solution 
of N chains of n blobs each: 



where 



= <S'i„tra(fc) + -2 X! X! ^"/3(^)' 
a (3 

N 



(3) 



^im^^pO^- 



i^r a=l 



(4) 



and fia is the position of the centre of blob a on chain 
i] Sai3ik) is the partial intermolecular structure factor 
for blobs a and /?, while S'intra(fc) is the intra-molecular 
structure factor of one polymer. Consider first an iso- 
lated polymer of n blobs {S{k) = 5'intra(fc))- Results 
for the multi-blob model with n = 62 and n = 602 are 
plotted versus fc • i?g on a log-log scale and compared to 
the structure factor of a full monomer representation of 
a SAW with 1000 monomers in figure [21 All the data 
follow the power law S{k) ~ a{kRg)~^/'^ in the scaling 
regime 1 < kRg < Rg/B, where B is the effective seg- 
ment length (rg in the multi-blob representation and b 
in the full monomer representation), with i' ~ 0.5988, 
close to the expected value of the Flory exponent, and 
a ~ 1.15, close to the renormalization group prediction 
14| . The multi-blob representation does preserve the cor- 
rect scaling behaviour of the single chain structure factor. 
The structure factor of the full monomer representation is 
approximately related to that of a single blob, Sbiob(fc), 
and to the structure factor of a chain of blobs, Sc{k) 
by the factorisation approximation, S{k) = Sb{k) ■ Sc{k) 
[iSj.The small k limit of this relation yields an estimate 
of the radius of gyration Rg of a SAW in terms of those of 
the chain of blobs {Rgc) and of a single blob (r^), namely 

d2 _ 7-)2 I 2 

We next consider S{k) in the semi-dilute regime. 
Structure factors obtained from MC simulations of 108 
and 500 chains of n = 5, 10, 20 blobs are plotted in figure 
3 for p/p* = 1.15 and 9.23 versus kRg, for a fixed overall 
length L of the polymer (i.e. a fixed value of Rg). The 
overall agreement between the three multi-blob represen- 
tations is seen to be good for krg < 1 {kRg < n") at the 
lower density, as one might expect within a consistent 
coarse-graininig procedure. The significant discrepancy 
observed at the higher density illustrates the fact that for 
p/p* — 9.23, the minimum number n of blobs needed to 
satisfy the requirement pb < pi is n — 18. Figure 4 shows 
the inter-molecular part of the total structure factor ([3]), 
as calculated with n = 5, 10 or 20 blobs; the agreement 




FIG. 2: (Color online) Intramolecular structure factor 
<S'intra(fc) of a singlc SAW polymer versus kRg on a log-log 
scale. Squares: MC data for an L = 1000 SAW polymer; 
closed circles and open circles: MC data for n = 62 and 
n = 602 blob representations. Full line: scaling regime 
1.15(fc_Rg)^'''^, with f ~ 0.5988: dashed line: Debye structure 
factor for a Gaussian chain. 
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FIG. 4: (Color online) Intermolecular structure factor from 
MC simulation of m = 5 (circles), 10 (squares) and 20 (tri- 
angles) blob representations of SAW polymers at densities 
p/p* = 1.15 (closed symbols) and 9.23 (open symbols) versus 
kRg. The dashed curves are predictions of the RPA derived 
from expression ([5]), using the intramolecular form factor (for 
n = 10). The RPA is seen to be accurate for kRg < 2.5, but 
to underestimate short-range correlations {kRg > 2.5). 
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FIG. 3: (Color online) Total (open symbols) and intramolec- 
ular (closed symbols) structure factors versus kRg from MC 
simulations of n = 5 (circles), n = 10 (squares) and n = 20 
(triangles) blob representations of SAW polymers at p/p* — 
1.15 (left frame) and 9.23 (right frame). The straight lines 
correspond to the scaling regime as in fig. [2] 



reads: 



S{k) = 



Lu{k)/n 
1 — pbuj{k)c{k) 



(5) 



Since the effective interaction between blobs is soft (cf. 
([T])), the direct correlation function is well approximated 
by the RPA closure c(fc) ~ —v{k)/kBT, where hats 
always refer to Fourier transforms. Gathering results: 



S{k) 



){k)/r 



1 + PbCj{k)v{k)/kBT 



(6) 



which generalises the classic RPA results for the structure 
of polymer solutions [3 12) to multi-blob chains. 
The advantage of the result ([5]) is that the RPA on which 
it relies is accurate because of the weakness of the blob- 
blob interaction. Using ujik = 0) = n and the Fourier 
transform of ([T]) one arrives at the compressibility: 



between the various coarse-grained representations is ex- 
cellent. The fc ^ limit of the structure factor yields 
the osmotic compressibility xt via pksTxT = S{k 0). 
In the same limit all partial intermolecular structure fac- 
tors become identical so that the compressibility is also 
given by pksTxT = 1 + s„/3(fc ^ 0). The PRISM ap- 
proximation assumes all partial structure factors to be 
identical for all /c, i.e. Sapik) — s(fc) ^ a(3, which is only 
true in the limit of very long polymers, where end effects 
may be neglected [12|. PRISM then relates the total 
structure factor to the form factor Lb{k) — nSintia.{k) of 
the polymer, and to the Ornstein-Zernike direct corre- 
lation function c{k). For chains of n blobs the relation 



pksTxT 



1 + pbuA (^) 



3/2 „3 ■ 

9 



(7) 



For given p and polymer length L, the number n of blobs 
is chosen such that pb = np — pi, i.e. n ~ (p/p*)^^'"'^" 
Hence: 



pksTxT 



1 



(8) 



where C = iAiirlaf'^ /At: ~ 3.23. For p » p* ® 
leads back to the des Cloizeaux scaling for the osmotic 
pressure H ~ {p/ p*f''/^^''-^'> - {p/p*f'^ [1]. Numerical 
values obtained for the compressibility from eq ([H]) agree 
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with MC data (as shown in figure 3) within a few per- 
cent, except at the highest density where a discrepancy 
of 20% is probably due to large statistical uncertainties 
in the MC data. Another consequence of the weakness of 
the effective interaction between blobs is that solutions 
of multi-blob chains lend themselves readily to thermo- 
dynamic perturbation theory {l9] . Choosing a system of 
non-interacting Gaussian coils of length n and spring con- 
stant k = SkBT/Bg as reference system, the free en- 
ergy _F of a solution of an equal number of n-blob chains 
satisfies the Gibbs-Bogoliubov inequality (191] 



F<Fo + (Vn~Vjl"^)^ 



(9) 



where Fq is the free energy of the reference system, Vn 
and Vn^'' are the total energies of the interacting n-blob 
system and non-interacting reference system, and the av- 
erage is taken with Gaussian statistics. The two terms on 
the right hand side are readily calculated analytically as 
functions of the bond length Bq of the Gaussian chains. 
The right hand side of ^ is then minimised with re- 
spect to Bq to provide the best variational estimate of 
F. The optimum Bq leads directly to an estimate of 
the radius of gyration of interacting polymers in a semi- 
dilute solution, namely Rg = n°-^_Bo/V6 = 
where rg ~ {L/nY is the radius of gyration of one blob 
and bo — (Bo/rg) is determined by the free energy min- 
imisation; bo increases with n and saturates at a value 
bo ~ 2.24 for n > 100. For a given polymer density p and 
length L, n is once more determined by the requirement 
that the blob density pb = np equals its overlap value pi ; 
this leads to the desired result: 



Rg = cbL" 



P 



(10) 



0.40 for sufficiently large n. Eq. (fTU)) shows 



correct slow decrease of Rg with density {Rg 



where c 
the 

(p/p*)~°-^^^) predicted by scaling theory [if. Note that 
in the melt limit n ^ L, Rg scales like L^/^, as expected. 

The present multi-blob coarse-grained description pro- 
vides a fresh look at semi-dilute polymer solutions and 
allows earlier work on dilute solutions P, [3, S 0, S IB] to 
be extended to concentrated solutions of long chains, us- 
ing transferable, state-independent pair interactions be- 
tween blobs. The softness of the latter, and reduction of 
the number of degrees of freedom, lead, on one hand, to 
a speed-up of simulations by orders of magnitude, thus 
allowing the simulation of large numbers of high molecu- 
lar mass polymers, while, on the other hand, opening up 
the systematic use of standard methods of statistical the- 
ories of fluids. We have shown that the coarse-graining 
procedure preserves standard scaling laws of semi-dilute 
polymer solutions, while providing reliable estimates of 



amplitudes (i.e. numerical pre- factors). We are presently 
extending the multi-blob methodology to more general 
homopolymer models involving effective attractions be- 
tween adjacent monomers (and hence solvent quality or 
temperature effects) , to brushes of grafted polymers [20| , 
to colloid-polymer mixtures and to the self-assembly of 
block copolymer solutions, going beyond the minimal 
two-blob representation of diblock copolymers used in 
references [ig. \vi\ . 
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